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Ordering of the classical Heisenberg antiferromagnet on the triangular lattice is studied by 
means of a mean-field calculation, a scaling argument and a Monte Carlo simulation, with 
special attention to its vortex degree of freedom. The model exhibits a thermodynamic transi- 
tion driven by the Z2-vortex binding-unbinding, at which various thermodynamic quantities 
exhibit an essential singularity. The low-temperature state is a "spin-gel" state with a long 
but finite spin correlation length where the ergodicity is broken topologically. Implications 
to recent experiments on triangular-lattice Heisenberg antiferromagnets are discussed. 
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Ordering of geometrically frustrated magnets has attracted much recent interest. Anti- 
ferromagnetic Heisenberg model on the two-dimensional (2D) triangular lattice is a typical 
example of such geometrically frustrated magnets. Interest in the ordering of the model has 
been enhanced by recent experiments on various triangular-lattice Heisenberg antiferromag- 
nets (AFs) including, 5=3/2 NaCrOa,^'^) S=l NiGa2S4^"^) and 5=1/2 organic compounds, 
K-(BEDT-TTF)2Cu2(CN)3^) and EtMe3Sb[Pd(dmit)2]2-^^ While most of these compounds ex- 
hibit a spin-liquid-like behavior at low temperatures without the conventional magnetic long- 
range order (LRO), all of them turn out to exhibit a weak but clear transition-like anomaly 
at a finite temperature. This "transition" seems intrinsically 2D, being neither the standard 
AF nor the spin-glass transition. Thus, it remains most interesting to clarify the nature of the 
ordering process of the frustrated AF Heisenberg model on the triangular lattice. 

It is now established that the Heisenberg AF on the triangular lattice interacting via 
the nearest-neighbor bilinear interaction exhibits a magnetic LRO at T = 0, the so-called 
120-degrees structure, in either case of quantum S = 1/2 spin^*^'^^-* or classical S = 00 spin. 
Because of the two-dimensionality of the lattice, the magnetic LRO is established only at 
T = 0, while the associated spin correlation length diverges exponentially toward T = 0. 
Some time ago, it was demonstrated by one of the present authors (H.K.) and Miyashita that 
the model bears a topologically stable point defect characterized by a two- valued topological 
quantum number, a Z2 vortex. Kawamura and Miyashita suggested that, in contrast to its 
unfrustrated counterpart, the triangular Heisenberg AF might exhibit a thermodynamic phase 
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transition at a finite temperature driven by the binding-unbinding of the Z2 vortices, where 
the vortex correlation length ^t,, corresponding to the mean separation of free Zi vortices, 
diverges keeping the spin correlation length ^ finite. 

It should be stressed that the proposed vortex transition is of different character from 
the Kosterlitz-Thouless (KT) transition of the 2D XY model in that there appear two dis- 
tinct length scales, or two different stiffnesses (energy scales), in the former: One is a vortex 
correlation length diverging at the vortex transition T = T^, and the other is a spinwave 
correlation length ^gyj staying finite at T = T^. Counterview to such a picture is that there is 
no such "vortex-spinwave decoupling" that the Z2-vortex transition is merely a rapid crossover 
rather than a true thermodynamic transition. 

Under such circumstances, we wish to investigate in the present Letter the nature of 
the Z2- vortex ordering of the classical AF Hcisenberg model on the triangular lattice, first 
analytically by use of mean-field and scaling analyses, and then numerically by use of a Monte 
Carlo (MC) simulation. On the basis of these analytical and numerical results, we shall discuss 
the recent experimental results on various triangular AFs. 

The model we consider is the classical Heisenberg AF on the two-dimensional triangular 
lattice, whose Hamiltonian is given by "H = J X^<jj> Si ■ Sj, where J > and the sum is taken 
over all nearest-neighbor pairs. 

As demonstrated in ref.,"'^^-' the Heisenberg spins, or more precisely the chirality vectors, 
circulate around a vortex core making a topologically stable vortex, whereas whether they 
circulate in clockwise or counter-clockwise fashion docs not make clear distinction topologi- 
cally. At T = 0, an isolated vortex of its radius R has an energy Ey ~ clogi? -|- /i, with c an 
energy constant and /i the vortex core energy. Stability of a single vortex against the cntropic 
effect yields a rough estimate of the vortex transition temperature = c/2 above which an 
isolated Z2 vortex is spontaneously generated. 

Suppose that there exist free vortices on the lattice of size N = Lx L. Mean separation 
between free vortices is given by = \j ^fn^ where n^, = N^/N is the number density of free 
vortices. Neglecting the correlation effects between vortices, the free energy of an assembly of 
vortices might be given by F„ ~ (cln,^t, -|- iJL)Ny — T\a[N\/{Ny\{N — Ny)\)]. By minimizing 
with respect to ra^, we get n„ ~ exp[— (^u — Ty)/(T — T^)] with Ty = c/2 as above. Then, the 
vortex correlation length ^„ corresponding to the mean separation of free vortices is obtained 
as « exp[2^^^^-j]. Toward T = T^, ^„ diverges exponentially, and free vortices disappear 
at T < r„. Previous MC suggested ~ 0.28 (in units of J)^^) and /x ^ 1.65.^2) 

In the above mean-field analysis, we have neglected the correlation effect (or the screening 
effect) between vortices. In reality, the screening effect would make neighboring vortices getting 
closer in distance. In view of this, we generalize the mean-field expression of to the form 

e.~exp[(A/(r-r,))"], {T>Ty), (1) 
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where A is a constant and the exponent a is expected to be smaller than the mean-field value 
unity, i.e., a < 1. While one has a = 1/2 in the standard KT transition, precise criticality 
of the Z2-vortex transition has not been elucidated yet. In any case, since the standard 2D 
Coulomb-gas description is not directly applicable to the Z2-vortex ordering, there seems to 
be no strong reason to believe that the exponent a is exactly a = 1/2. 

From a general scaling argument, the singular part of the free energy around is given by 
fv ~ ^v'^- At T = T^,, /t, exhibits an essential singularity. Though weak, an essential singularity 
is a well-defined singularity so that the expected Z2-vortex transition is a thermodynamic 
transition, not just a rapid crossover nor a smeared transition. The specific heat is given by 

C ^ {A/{T - T„))2("+i) exp[-2 {A/{T - T„))"] 

+ [regular part], (T > T^), (2) 

where the regular part represents the contribution from degrees of freedom other than free 
vortices such as spinwaves. While the specific heat exhibits only a weak essential singularity at 
T = T„ where the first vortex-pair unbinds, it would exhibit a non-singular peak or a shoulder 
at a temperature slightly above where the vortex-pair unbinding occurs most extensively. 
The precise form of the specific heat, however, would largely depend on the form of the regular 
part, which is a non-universal property sensitive to the details of each particular system. 

Spin correlation function C(rjj) =< Si ■ Sj > is affected both by vortices and spinwaves. 
We assume here that the spin correlation can be factorized into the vortex part Cu(r) and the 
spinwave part Cs„,(r), C(r) ~ Cv{r)Csw{'>^), as in the KT ordering of the 2D XY model. 
Assuming the standard exponential decay for each correlation function at long distances, 
C(r) exp(-|), Cy{r) exp(-^) and Csw{r) « exp(-^), we get ^ = CvCsw/iCv + Csw)- 
Note that the vortex correlation length ^t,, diverges at T = T^, , whereas the spinwave correlation 
length ^sw stays finite at T = Ty, which monotonically increases with decreasing T and 
eventually diverges toward T = 0. In the high-temperature regime T » where « ^sw, 
spin correlation is dominated by vortices, while in the low-temperature regime T « T„, 
dominantly by spinwaves. At temperatures slightly above where >> ^sw, C is given by 

with given by Eq.(l). Hence, ^ exhibits a weak essential singularity at T = T^, remaining 
finite there. The relative importance of vortex and spinwave excitations is interchanged around 
a crossover temperature Tx > at which = ^sw The appearance of two distinct length 
scales and ^sw, or two distinct coupling constants (stiffnesses), is a unique feature of the 
Z2-vortex ordering of the Heisenberg system. This is in contrast to the standard KT ordering 
of the XY system where there exist only one length scale (one stiffness). 

The low-temperature phase realized below is an "ordered" state in the sense that 
it cannot be reached adiabatically from the high-temperature phase, though it is a "spin 
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paramagnetic" state with a finite spin correlation length. It is essentially a spinwave state 
where the full ergodicity is broken topologically since only the vortex-free sector is allowed on 
long length scale in the phase space. Dominance of spinwave excitations below would give 
rise to the specific heat proportional to at T < T^,, once the quantum effect has been taken 
into account (in the pure classical model, the specific heat is not realized). ""^^^ Although 
the spin correlation length ^ or the spin correlation time r still remains finite even below Ty, 
the suppression of free vortices might lead to a substantial growth of ^ or r at and below Ty, 
leading to a near-critical (but not true critical) state. We call such a topologically ordered 
state with a finite but long spin correlation length and correlation time a "spin gel" state. 
This state differs in nature either from the standard AF ordered state, from the spin-liquid 
state, or from the the spin-glass state. 

Finiteness of the spin correlation length ^ at T = T^, has important consequences on the 
response of the system against weak perturbative interactions. For example, a magnetic field, 
which couples to the spin via the Zeeman term, reduces the Hamiltonian symmetry so that 
one may suspect that even an infinitesimal magnetic field suppresses the Z2-vortex transition. 
This is not necessarily the case, however, in the situation where ^ stays finite at T = Ty. 
If H is sufficiently small satisfying HS^{T = T^,)^ ^ ksTy, the Z2-vortex transition and the 
low-temperature "spin gel" state might remain essentially the same as in the zero- field ones. 

Similar situation is expected for other perturbative interactions such as the magnetic 
anisotropy D or the interplanar 3D coupling J'. If these perturbative interactions are weak 
enough satisfying ^ ksTy or J'^^ S ksTy, the Z2-vortcx transition and the low- 

temperature spin-gel state might well persist. Tims, if ^ stays sufficiently short at T = Ty, the 
Z2-vortcx transition and the spin-gel state of the 2D Hcisenberg system might be "protected" 
against the weak anisotropy and/or the interplanar interaction. 

Under an infinitesimal field H, the singular part of the free energy still keeps the same 
form as the zero-field one where the transition temperature Ty and the constant A now weakly 
depend on H but in a regular way, Ty{H) = T^{0) + cH^ + ■■■ and A{H) = A{0) + c'H^ + ■■■. 
Prom this, prediction for the zero-field susceptibility x follows, 

X ^ (A/iT - Ty)r+' exp [-2 (A/(T - T,))"] 

+ [regular part], (T > T^), (4) 

where the exponent a and the constant A are common with the ones in Eqs.(l) and (2). Again, 
X exhibits a weak essential singularity at T = T^. Similar expression can also be derived for 
the nonlinear susceptibility xs , where the exponent in the prefactor of Eq. (4) 1 + a should be 
replaced by 2(a -|- 1). 

In order to examine the validity of the above analysis, we also perform a MC simulation 
of the model based on the standard heat-bath method combined with the over-relaxation 
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Fig. 1. (Color online) Temperature dependence of the specific heat (a) and the susceptibility (b) per 
spin. The curves are the fits based on Eqs.(2) and (4). The inset is the Curie- Weiss plot for the 
inverse susceptibility. 



The specific heat C and the susceptibility x ^-I's shown in Fig.l. The specific heat exhibits 
a distinct but rounded peak around T ~ 0.318, located shghtly above T„. More precise T^- 
value will be determined below in Fig. 3 as ~ 0.285 it 0.005, together with the a-value, 
a = 0.42 lb 0.15. The Curie- Weiss temperature Tew is estimated as Tew — 2.32: See the inset 
of Fig. 1(b). The C data are then fitted to our theoretical expression Eq.(2) where the regular 
part is approximated by the fourth-order polynomial of T. As can be seen from Fig. 1(a), the fit 
works well with A = 0.66±0.05. The x data can also be fittable to Eq.(4) with A = 0.74±0.55, 
which agrees with the A value determined above from C. 

In Fig. 2, we show on a semi-log plot the finite-size spin correlation length associated 
with the 120-degrees structure. A scaling relation ^2l/Cl ^ f{£,L/L) is expected to hold 
for larger L}^^ As shown in the inset, the data turn out to scale well in this form for larger 
L, yielding a scaling function /(x). Then, following ref.,^®) we extrapolate to the bulk 
correlation length ^ = S^^o on the basis of the above scaling relation, the resulting ^ being 
given in Fig. 2. The change in the behavior of ^(T) is discernible around T = T^. The data 
can be well fitted to our theoretical expression Eq.(3) with A ~ 0.94 it 0.4 and a = 0.42 
(fixed), where we fit In ^sw by the fourth-order polynomial of T. Note that, in contrast to 
the ^ data reported in ref.,^^^ our ^ data, or its temperature derivative, does not show any 
appreciable discontinuity at T = T^, exhibiting only a weak singularity consistent with an 
essential singularity. This difference from ref.^^^ probably comes from the fact that the two- 
loop RG formula used in ref.^^-* in performing an L = oo extrapolation is no longer valid 
around T = T^. 

In the present model, although the spin correlation length remains finite at T = T^, it 
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Fig. 2. (Color online) Temperature dependence of the finite-size spin correlation length, together 
with the bulk one. The inset is a scaling plot, i.e., ^2l/Cl plotted vs. (,lIL. 



becomes quite large, ^ ~ 1700. Our fitting analysis also gives an estimate of the crossover 
temperature Tx ~ 0.294 and the crossover length /x — 380. It should be noticed, however, 
that the value of ^ at T = (at T = Tx ) is a non-universal property depending on the details 
of each system. 

In the previous studies, as an order parameter characterizing the Z2-vortex transition, 
either a Wilson- loop (vorticity function) or a vorticity modulus^*^^ has been proposed. Here, 
we calculate following Refs."'^^'^"'^^ the latter quantity, the vorticity modulus w, defined as the 
free-energy cost against a vortex formation AF divided by InL, v = AF/lnL. In Fig. 3, we 
show the temperature derivative of the vorticity modulus, -dv/dT, calculated from appropriate 
fluctuations. The data exhibit a sharp peak which sharpens with L, suggesting the occurrence 
of a phase transition. In the inset, the peak temperature Tp(.ak{L) is plotted versus l/ln(L). 
Indeed, Eq.(l) implies a relation Tpeafcl^) -7^ ~ (l/lnL)^/". This fit yields = 0.285±0.005 
and a = 0.42 it 0.15, which are the and a values quoted above. 

We now wish to discuss possible implications of our results to recent experiments on several 
triangular-lattice Heisenberg antiferromagnets. We note 5=3/2 NaCr02 {Tqw — 290K)^'^) 
and S=l NiGa2S4 {Tew ~ BOK)^"'') exhibit a strikingly similar ordering behavior in spite of 
the difference in their integer /half- integer spin quantum numbers. Indeed, NaCr02 (NiGa2S4) 
exhibit a clear but rounded specific-heat peak at Tpgafc = 41K (lOK), whereas a transition- 
like anomaly is observed at Tf ~ 30K (8.5K), a temperature slightly below Tp^ak, where the 
spin dynamics is rapidly slowed down giving rise to a quasistatic internal field. However, the 
spin dynamics probed by NMR, NQR, ESR and /xSR is not completely frozen, but the spins 
remain slowly fluctuating even below Tf, unlike the conventional AF or the spin glass. Such 
a dynamically fluctuating ordered state extends over a wide temperature range, down to lOK 
(2K). The spin correlation length determined from neutron scattering is kept finite even at 
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0.35 

Fig. 3. (Color online) Temperature dependence of the minus of the temperature derivative of the 
vorticity modulus. The inset represents the size dependence of the peak temperature. 

and below Tj, i.e., ~ 20 (5) lattice spacings. At least in NiGa2S4, dynamical freezing at 
T = Tf accompanies a weak structure in the susceptibility, while application of magnetic 
fields greater than 100 ~ 1000 G gradually changes the slowly fluctuating state into more 
conventional frozen state. 

These experimental features are fully consistent, at least qualitatively, with the Z2-vortex 
order, if the experimental freezing temperature Tf and the low-temperature state are identi- 
fied as and the "spin-gel" state. In particular, energetics seems appropriate: Our simulation 
yields Tpeafc/^CH^ — 0.137 and T^/Tcw ~ 0.123, which are close to the corresponding experi- 
mental values 0.14 and 0.10 (NaCr02), and 0.13 and 0.11 (NiGa2S4). Furthermore, in case of 
NiGa2S4, a weak anomaly observed in the specific heat and the linear and nonlinear suscep- 
tibilities are well fittable to our theoretical formula (2) and (4),^^^ while the crossover- field, 
~ lOOOG, estimated from our theoretical formula with the experimental value of ^ ~ 8 seems 
consistent with the //SR measurement. Major quantitative discrepancy between experiments 
and our present simulation lies in the magnitude of the spin correlation length at the transi- 
tion. In our simulation, ^ is of order 1000 near T = T^, while it is only 5 or 20 experimentally. 
So, some mechanism, not taken into account in the simplest classical Heisenberg model, is 
required to explain the shortness of As emphasized, the magnitude of ^ is a nonuniversal 
property governed by the non-vortex physics, which could significantly be reduced, say, by 
quantum fluctuations, charge fluctuations, or further frustration effects associated with the 
interaction other than the main exchange coupling. 

S'=l/2 organic AFs k-(BEDT-TTF)2Cu2(CN)3 and EtMe3Sb[Pd(dmit)2]2 also exhibit 
a spin-liquid behavior without the conventional magnetic LRO down to low temperature. 
Recent measurements revealed a weak anomaly at a finite temperature, 4~ 6K for BEDT- 
TTF^) and ~ 5K for dmit,^) where the NMR Ti exhibits an anomaly (or the spectra broaden) 
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and the specific heat exhibits a broad hump. Wc note that the Z2- vortex order is one possible 
candidate of the experimentally observed anomaly. Of course, for the vortex to be a meaningful 
excitation, minimal amount of the noncollinear spin short-range order of a few lattice spacings 
is required, which, however, is not unlikely in these compounds. In order to substantiate this 
conjecture, further study is required. 
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